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Abstract. We study the entanglement dynamics in the externally-driven single- mode Dicke 
model in the thermodynamic limit, when the field is in resonance with the atoms. We compute 
the correlations in the atoms-field ground state by means of the density operator that represents 
the pure state of the universe and the reduced density operator for the atoms, which results from 
taking the partial trace over the field coordinates. As a measure of bipartite entanglement, we 
calculate the linear entropy, from which we analyze the entanglement dynamics. In particular, 
we found a strong relation between the stability of the dynamical parameters and the reported 
entanglement. 
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1. Introduction 

Understanding the properties displayed by the entanglement of physical systems is one of the 
fundamental purposes of quantum information theory [T]. A quantum system composed of two 
or more entangled subsystems has the interesting property that although the state of the total 
system can be well defined, it is impossible to identify individual properties for each one of 
its parts. A subject that has attracted recent interest is linked to the relationship between 
entanglement and certain physical properties of many body quantum systems; this has taken 
special interest in relation to quantum phase transitions (QPTs) (2j El 0]. A QPT in a many 
body system strongly influences the behaviour of the system near to the critical point, with 
the consequent appearance of long-range correlations in the ground state. Although several 
proposals exist, currently there is no complete theory of multipartite entanglement, and the 
common techniques are based on bipartite decompositions of the total system. This type of 
decomposition has made it possible to study the entanglement in the Dicke model in thermal 
equilibrium j5j[6l[8]. Here, we study the entanglement properties of a many-body system in the 
non-equilibrium regime, in which the total system has a unitary dynamics but is not an isolated 
system. For this purpose we study the single mode "externally-driven" Dicke Hamiltonian 



where J z = Yli=i J* j J± = Yli=i J± are collective atomic operators, loq is the level splitting of 
the atoms, uj is the frecuency of the bosonic mode and g{t) = g + Ag cos fit is the time dependent 





atom-field coupling. We assume that Ag is a fraction of the static coupling g. An exact 
diagonalisation of the problem Eq. (pQ) has been previously carried out by means of the Holstein- 
Primakoff transformation for the case of a static atom-field coupling g(t) = g = constant [5j[6], 
but there is no known solution for the time dependent atom-field coupling case. 



2. The reduced density operator 

When studying a composed quantum system whose dynamics is unitary, it is interesting to study 
the dynamics of the parts of the system. In contrast to the dynamics of the whole system, the 
dynamics of the subsystems is not unitary. However, it is possible to make a description of 
the subsystem through the reduced density operator. The Schrodinger equation for the single 
mode "externally-driven" Dicke Hamiltonian has an exact solution when the field is in resonance 
with the atoms co = loq. With the purpose of formulating the exact solution, we introduce an 
abstract coordinate representation through the coordinates of the physical field x and the atoms 
coordinate y. In terms of these coordinates, the ground state of the universe (atoms+field) is 
given by the expression [9] 
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It is interesting to note that the dynamics of the system is influenced by the dynamics of the 
auxiliary parameters B^(t), which are solutions of the Mathieu equation [10] 



B^(t) + [(1 =F 2g) =F (2Ag) cos nt]B^(t) = 0, (3) 



subject to the Wronskian condition 

B^(t)(B^)*(t) - B*(t)(B*)*(t) = i. 

We describe the dynamics of the total system in terms of the density operator, where the 
pure state of the universe is represented by the operator pc(t) = |^o_,o +5 *)(^'o_,o +5 *|- I n the 
coordinate representation, the density matrix takes the form 



PG (x',y';x,y,t) = %_ fi+ (x,y,t)^ _ fi+ (x',y',t). (4) 

We establish a bipartite decomposition of the universe and study the reduced dynamics of the 
atoms through their reduced density matrix (RDM). In so doing, we calculate the partial trace 
over the physical field coordinate x, 
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Figure 1. Phase space trajectories of the auxiliary dynamical parameters B^(t) for Ag = O.lg, 
and a value of the static coupling in a) the stable region g = 0.40, b) the stable region g = 0.46, 
and c) the unstable region g = 0.38. 



The ground state of the universe *o_,o+ {x, y, t) given in Eq. ([2]) is Gaussian, which facilitates 
the integration of Eq. ([5]). After some algebraic calculations, we obtain the following result for 
the reduced density matrix 
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where c = s = l/\/2 and 3?e(^ =F ) is the real part of the function ^(t) = — . In the 
particular case Ag = (g(t) = g), and uj = ujq, this function is independent of time and 
becomes £ T (t) = e T = ^(1 =p 2p). Hence, our model reproduces the result previously reported 
in the literature for the atom's reduced density matrix for the single mode Dicke model in 
thermal equilibrium [5j [6] . 



3. The linear entropy 

We consider a pure state of a composed system (the universe), say system AB, represented by 
the density operator pab- The linear entropy for the reduced density operator pa = trg (pab) 
of subsystem A is defined by La = 1 — tr^p 2 ^), where trA (trs) denotes the trace over the 
subsystem A (B) [7J. This gives a measure of purity of the reduced density operator pa (o pb) 
of one part of the total system (in the bipartite decomposition of the universe). If the pure state 
of the universe is separable, the reduced density operator of one part of the system represents 
a pure state and, as a result, its linear entropy must be zero. Similarly, if the pure state of the 
universe is a maximally entangled state PQ, the linear entropy equals 1/2. In the context of our 
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Figure 2. Linear entropy as a function of time for periodic atom-field coupling with Ag = O.lg; 
a) g = 0.40, b) g = 0.46, and c) g = 0.38. 



externally driven single mode Dicke model, it is possible to use the linear entropy as a measure 
of bipartite entanglement for the atoms- field system. The linear entropy for the atoms' time 



dependent reduced density operator, Eq. ([6]), is given by L(t) = 1 — tr 
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. By using 



the density matrix representation of Eq. §§§ , we explicitly calculate the linear entropy as 
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which gives, after some algebraic calculations, the result 

L(t,g) = l 
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where the parameters a, (5 and A are defined as 
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4. Results 

Interestingly, the linear entropy exhibits a time dependence that is not determined by the global 
phase of the ground state wave function of the universe, Eq. ([2]), but by the auxiliary dynamical 
parameters B^(t). In the particular case Ag = O.lg, and Q = 1, the stability zones of Eq. ([3]) 
are known [9], and the study of the corresponding solutions is based on the Floquet's theorem 
for second order differential equations with time periodic coefficients. For this, the solutions 
of Eq. © have the general form B ± {t) = exp(iF ± t)4> ± (t), where ^(t + T) = ^(t) and F ± 
is the Floquet exponent, which depends on the shape of the driving. For driving functions for 
which is complex, the solution becomes unstable. In the stable regime, F^ is a real number 
[TO] . In Figs. 1(a), (b) we consider the phase space representation of the trajectories of the 
auxiliary dynamical parameters (t) . For values of g that belong to a common stability zone 
(g = 0.4 and g = 0.46) [9] the solutions are bounded in the phase space. In order to establish 
a relationship between the stability and the entanglement dynamics, we consider the linear 



entropy in Figs. 2(a), and (b) for the values of the static coupling described above. Our choice 
of the initial conditions for Eq. ([3]) implies that the system starts in thermal equilibrium, but 
is described by a non-separable quantum state. Strikingly, the dynamics exhibits a behaviour 
whereby the system experiences a disentanglement process, before reaching a maximum value 
of entanglement, and then successive collapses and revivals. 

Figure 1(c) shows another interesting behaviour: when the static coupling belongs to the 
instability zone {g = 0.38) of B~(t), the solution B + {t) is bounded while the solution B~ is 
unbounded in the phase space. For this value of the static coupling, the linear entropy oscillates, 
before it reaches the stationary state with linear entropy L = 1, as shown in Fig. 2(c). In [9j 
we discuss the relationship between the stability zones and the localization of the ground state 
wave function of the universe. In order to study this relation, we define the characteristic length 
P(t) = v2|-B-F (Oli which is bounded in the stable zones and unbounded in the unstable zones. 
The numerical simulations of the ground state probability density show that for values of the 
static coupling that belong to the common stability zones the probability is localized in the 
abstract x — y space and presents an oscillatory behaviour. In contrast, for values of the static 
coupling that belong to the unstable zones, the probability density has an oscillatory behaviour. 
However, when the system evolves, the probability density is systematically stretched in a fixed 
direction, with a consequent dilation in the perpendicular direction. 

5. Conclusions 

We have obtained exact results for the atoms' reduced density operator and the linear entropy 
for the externally driven Dicke model, in the thermodynamic limit. These results allow us to 
study the entanglement dynamics through the auxiliary dynamical parameters. In contrast, 
in the Dicke model in thermal equilibrium (5J [6], the atoms-field entanglement does not have 
dynamics because the temporal dependence of the ground state of the universe is given by a 
global phase of the wave function. In this work, the auxiliary dynamical parameters and its 
stability properties determine the entanglement dynamics of the system and the behaviour of 
the universe ground state wave packet. 
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